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Abstract 

The study of latitudinal body size clines can illuminate processes of local adaptation, but there is a need for an increased 
understanding of the relative roles of genetic variation, environmental effects, and interactions between them, in shaping 
such geographical patterns. For this reason, we combined an investigation of a museum collection of the common blue 
butterfly Polyommatus icarus (Rottemburg) (Lycaenidae: Polyommatini) from Sweden with a common-garden experiment 
in the laboratory, using strains reared from individuals collected from three different latitudes. Sizes of the field-collected 
butterflies tended to smoothly decrease northwards in a latitudinal cline, but suddenly increase at the latitude where the 
life cycle changes from two to one generations per year, hence allowing more time for this single generation. Further 
north, the size of the field-collected butterflies again decreased with latitude (with the exception of the northernmost 
collection sites). This is in accordance with the "converse Bergmann" pattern and with the "saw-tooth model" suggesting 
that insect size is shaped by season length and number of generations along latitudinal transects. In contrast, under 
laboratory conditions with a constant long day-length there was a different pattern, with the butterflies pupating at a 
higher mass when individuals originated from southern populations under time stress to achieve a second generation. 
This is indirect evidence for field patterns being shaped by end-of-season cues cutting development short, and also 
suggests counter-gradient variation, as butterflies from the time-stressed populations over-compensated for decreasing 
larval development time by increasing their growth rates, thus obtaining higher mass. Hence, we found support for both 
adaptive phenotypic plasticity and local genetic adaptation, with gene-environment interactions explaining the observed 
field patterns. 
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Introduction 

The study of latitudinal life history trends within species 
has a strong potential to illuminate processes of local ad- 
aptation, which in turn can provide information on how 
populations may respond to factors such as climate 
change. It might be expected that geographical gradients 
in climate, and in length of the favourable season, should 
translate to biological gradients (for instance clines in an- 
imal life history traits) in a fairly predictable manner. To 
the extent that this is not the case would be evidence of 
our imperfect understanding of the processes behind loc- 
al adaptation. As long as we cannot predict trends in the 
characteristics of organisms along such relatively simple 
gradients, we have little hope of predicting or even un- 
derstanding biological variation in cases where the un- 
derlying environmental causes are less evident. That said, 
the study of latitudinal and altitudinal variation has 
proven to be a much more complex and challenging sub- 
ject than perhaps was originally thought when various 
"rules of thumb" regarding geographical variation was 
first proposed. 

One of the best known intra-specific geographical pat- 
terns in animals is "Bergmann's Rule", depicting a size 
trend with larger size at higher latitudes as is found espe- 
cially in the endothermic mammals and birds (Ashton et 
al. 2000). The underlying selection is in this case gener- 
ally assumed to be related to thermoregulation; the allo- 
metric relationship between body mass and surface area 
selects for bigger animals in cold areas. It has been sug- 
gested that larger size at higher latitudes in ectotherms 
may be due to non-adaptive responses to low temperat- 
ures (Van Voorhies 1996), but adaptive explanations can 
also be found to explain increased size with latitude, oth- 
er than thermoregulation. For instance, Boyce (1979) ar- 
gued that an increased degree of seasonality with latitude 
can select for increased body size. 

In insects and other arthropods the opposite size trend, 
i.e. decreasing size with increasing latitude (or altitude), 
has often been reported; this is the so-called "converse 
Bergmann's rule" (Mousseau 1997), that has been well 
documented (Mousseau and Rolf 1 989; Nylin and Svard 
1991; Blanckenhorn and Fairbairn 1995; Telfer and 
Hassall 1999; Johansson 2003). Still, some reports de- 
scribed an increase in size with increasing latitude also 
occurring in arthropods (Chown and Gaston 1999; 
Blanckenhorn and Demont 2004). There is clearly a 
need to explore the causes of this variation (Chown and 
Gaston 1999; Johansson 2003; Blanckenhorn and De- 
mont 2004). 

A possible explanation for the occurrences of variable 
patterns of body size with latitude in insects was sugges- 
ted by Chown and Gaston (1999). They emphasised that 
species with short life spans, inhabiting ephemeral 
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habitats, could be expected to increase their mass at 
higher latitudes. For such species the length of the fa- 
vourable season does not set the limits for development 
time, but instead e.g. the availability of food resources. 
Larger size in colder areas could thus follow. In contrast, 
insects with a generation length more similar in mag- 
nitude to the length of the season can instead be expected 
to be constrained by season length and decrease their 
mass at high latitudes and altitudes. These predictions 
were later supported in a review of Blanckenhorn and 
Demont (2004). 

The most elaborate theory of geographical size patterns 
in insects is that developed by Rolf (1980, 1983), and it is 
compatible with the hypotheses of Chown and Gaston 
(1999) and Blanckenhorn and Demont (2004). Rolf sug- 
gested that populations with many generations per year 
may often show an r-maximizing strategy (as is often im- 
plicitly assumed in insect optimality models), but that 
species or populations with few generations are better de- 
scribed as maximizing R y , the yearly rate of increase. 
The reason is that insects generally can spend the winter 
only in a specific developmental stage, so that whole gen- 
erations have to fit into the favourable season. Hence, 
Rolf predicted that when the available time for growth 
decreases, the insect will be selected to mature earlier at 
the cost of a reduced adult size that could also reduce 
fecundity. This scenario results in a detectable pattern in 
development time and size along an environmental 
gradient; with "either a monotonic increase with 'season 
length' ... or a 'saw-tooth' pattern (Roff 1983). 

The saw-tooth pattern (Figure 1) is predicted to occur 
when a species shifts strategy from, for example, bivoltine 
(two generations per year) to univoltine (a single genera- 
tion) within a study area. Insects in the bivoltine popula- 
tion near the transition zone to the univoltine zone live 
under time stressed conditions since their favourable sea- 
son has effectively been cut in two parts, with the theoret- 
ical result of short development time and small individu- 
als as elements of local life history adaptation. Con- 
versely, the univoltine populations near the transition 
zone to the bivoltine area will have surplus time, and 
hence can develop for a longer time period to reach a 
bigger size. Size patterns wholly or pardy consistent with 
such saw-tooth patterns have been described in some in- 
sect species (Masaki 1978; Mousseau and Roff 1989; 
Nylin and Svard 1991; Mousseau 2000; Johansson 2003). 
See also Petersen (1947; e.g. pp. 438-439) for what is 
probably the first mention of an idea fully in line with 
Roff s hypothesis, introduced to explain the complex lat- 
itudinal size patterns in the butterfly Pieris napi in 
Sweden. 

Observed field patterns can be due to complex interac- 
tions between Bergmann and converse Bergmann trends 
(Blanckenhorn and Demont 2004), further complicated 
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Figure I. Theoretical "saw-tooth" pattern. After Roff (1980, 1983) and Nylin and Svard (1991). 



by shifts in voltinism. In addition, there may be counter- 
gradient variation in growth rates (Conover and Schultz 
1995) which has the potential to erase field trends in size 
even though, or because, genetic differences are evident 
in common-garden experiments. Over- or undercom- 
pensating counter-gradient variation is a further possibil- 
ity and can result in field size trends that are difficult to 
predict (Blanckenhorn and Demont 2004). 

Butterflies are known to frequendy adjust their develop- 
mental time and growth rate according to photoperiod 
and temperature, and mass/size may also be a plastic 
trait (e.g. Wiklund et al. 1991; Leimar 1996; Gotthard 
1998; Fischer and Fiedler 2000). Plastic growth rates 
have the power to uncouple development time from final 
size (Wiklund et al. 1991; Nylin 1992; Abrams et al. 
1996; Davidowitz and Nijhout 2004) creating a three- 
dimensional time-rate-size relationship rather than the 
positive correlation between time and size explicidy as- 
sumed by Roff (1980, 1983). It is thus of special interest 
to determine whether insects with plastic growth rates, 
such as butterflies, still follow the patterns predicted by 
Roff and, if so, how such patterns are shaped by genetic 
or plastic variation among sites, or by gene-environment 
interactions. To this end, and as a test of the prediction 
that latitudinal trends in insects with generation lengths 
comparable to season length should show size trends ex- 
plained by season length (Chown and Gaston 1999), we 
examined collections of field-caught individuals of the 
common blue butterfly, Polyommatus warns (Rottemburg) 
(Lycaenidae: Polyommatini), and then reared strain from 
three different latitudes in a common environment. We 



found that in the field there is a pattern consistent with 
the saw-tooth model, but in the laboratory (with constant 
daylength and temperature) we found a different pattern, 
presumably resulting from counter-gradient genetic vari- 
ation. This suggests that field size patterns in this butter- 
fly are shaped by an interaction between genetic life his- 
tory differences and environmental end-of-season cues. 

Materials and methods 

Studied species 

The larvae of P. icarus feed on small herbs in Fabaceae, 
and the species spends the winter in a late larval instar. It 
is commonly found in grasslands and other open areas 
over the whole of the European continent, except for 
some islands on the northern tip of Scandinavia at latit- 
udes > 70° N (Dal 1978). Lycaenids are small butterflies 
generally considered to be relatively sedentary (Hughes 
2000). Therefore we assumed that P. icarus would have 
the opportunity to be locally adapted and show local ge- 
netic life history adaptations to climatic gradients, as was 
earlier found in some other lycaenids (Nylin and Svard 
1991). 

P. icarus is known to have two generations per year in 
Southern Sweden, e.g. on the Baltic island of Oland (loc- 
ality L in Figure 2; Professor O. Leimar, Stockholm 
University, personal communication). We initially as- 
sumed that the species is univoltine in Northern Sweden, 
with a transition zone in Central Sweden, but the details 
regarding voltinism were revealed in a study of the mu- 
seum collections of P. icarus (see below). P. icarus has been 
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Figure 2. Map of Sweden depicting the source populations for the common-garden experiment: L = low latitude 
(Oland), I = intermediate latitude (Stockholm), H = high latitude (Borgsjo). One longitude (degrees East, vertical line) 
and three latitudes (degrees North, horizontal lines) given. 



shown to respond to photoperiodic cues by plastically ad- 
justing its size and development time according to per- 
ceived date in the season (Leimar 1996). 



Museum collection 

The Swedish Natural History Museum has a collection of 
some 500 individuals of P. icarus collected from different 
latitudes in Sweden. The collection contains butterflies 
sampled between the years 1880 to 2001. Latitude and 
date of collection were determined from labels. Wing 
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lengths (base to tip) of all individuals were measured un- 
der a standard laboratory microscope (6x). The average 
of right and left wings was used. 

For purposes of illustrations and statistical analyses of size 
patterns, latitudes for collection sites were rounded off to 
nearest integer degree, whereas for purposes of studying 
the shift in generation number they were rounded off to 
nearest half degree. Rounding was done to increase 
sample sizes in latitudinal categories. We excluded the 
extreme latitudes of 55° and 69° due to very low sample 
sizes. 

Common-garden experiment 

Females were collected in the wild from three popula- 
tions for rearing of their offspring in the laboratory 
(Figure 2). Sixteen females were collected from the low- 
latitude population (two sites on Oland, 56.5° N), four 
from the intermediate population (Stockholm, 59.5° N) 
and two from the high-latitude population (Borgsjo, 
62,5° N). In Britain P. icarus shows a shift from a predom- 
inandy bivoltine to a univoltine phenology approximately 
at latitude 54°, but with variation among years (Asher et 
al. 2001). In Sweden the transition zone seems to be even 
more ambiguous and individuals following a bivoltine 
pathway can regularly be found at higher latitudes than 
in Britain (Professor O. Leimar, personal communica- 
tion, and see Results from the museum collection). 
However, the butterfly is clearly normally univoltine at 
the high-latitude site and generally bivoltine at the low- 
latitude site. In the intermediate area a second genera- 
tion is not uncommon but neither is it a dominant 
strategy. 

Length of the favourable season follows a strong and 
clear north-south pattern within Sweden (Nylin and 
Svard 1991) due to a combination of factors: the effect of 
latitude itself combines with generally higher altitudes 
and more inland climates in the north. The low-latitude 
and the intermediate sites are situated close to the Baltic 
Sea. The high latitude sampling area is more continental, 
about 80 km inland and for instance has a visibly differ- 
ent flora compared to the coastal regions at the same lat- 
itude (Nygren, personal observation). This implies that 
despite the uniform latitudinal distances between the tree 
sampling areas, the high latitude sample is actually cli- 
matically more distant to the intermediate area than sug- 
gested by the geographical distance alone. The location 
of the low-latitude sites on an island should even further 
enhance the climatic distances between sites, in terms of 
day-degrees available in the favourable season. 

To control for phenological changes in host plant quality, 
collected butterflies from all three sites were reared simul- 
taneously. Wild-caught females and later eggs from low 
latitudes were stored for a short time period in a cold 
room to delay larval hatching and achieve 



synchronization. The experiment was initiated over a 
ten-day period in the second half of July. Newly hatched 
larvae were individually placed in plastic jars containing 
fresh bird's foot trefoil, Lotus comiculatus L. (Fabales: 
Fabaceae) and reared in a climate chamber at a pho- 
toperiod of 22: 2 LD at 22°C (low latitude n = 320, inter- 
mediate n = 52 and high n = 34). The plants were 
checked daily and changed when signs of wilting were 
found. Larvae of different latitudinal origin were mixed 
in the climate room to prevent position effects, and also 
moved in the room at least every second day. 

Date of pupation was noted, and two days after pupation 
the pupal mass was measured and individuals were 
sexed. Wing length of adults was not taken for the 
laboratory-reared individuals, but pupal mass typically 
correlates well with adult weight and adult wing length 
among individuals within sexes in butterflies (Nylin 1992; 
Kemp 2000; Fric and Konvicka 2002) and moreover 
only differences among categories that were large enough 
to make the body measure chosen more or less irrelevant 
were used. 

Growth rate was not measured independendy of pupal 
mass, but for illustrative purposes a measure of growth 
rate was calculated as the natural logarithm of pupal 
mass divided by larval development time (in days). 

All statistical analyses were performed with 
STATISTICA 1999 Edition, Kernel release 5.5. 

Results 

Museum collection 

Above latitude 60° N, few P. icarus butterflies from a po- 
tential second generation, i.e. with a late collection date, 
could be found in the museum collection. At these latit- 
udes, collecting dates were distributed approximately 
normally with a peak at mid-July. At latitudes below 60° 
N, the distributions of collecting dates had two peaks (late 
June and early August). Consequently, the butterfly ap- 
pears to have shifted from a univoltine to a bivoltine life 
cycle at this latitude in Sweden. However, below latitude 
60° the butterfly probably is only partly bivoltine as not 
all individuals develop direcdy, and this generation may 
not occur in all years at all latitudes (the collection was 
not large enough to permit further analysis of phenolo- 
gical details). It is unclear where an obligate bivoltine life- 
style takes over, if ever within Sweden, but below lat. 
57.5° the second peak was consistendy large, and hence 
bivoltinism appears to be the dominant strategy. These 
results confirmed our initial assumptions regarding where 
the bivoltine and univoltine populations would be found 
in the field. 

In the "bivoltine area", below latitude 60° N, the butter- 
flies were classified as belonging to either the first or the 
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Figure 3. Wing length of Polyommatus icarus collected at different latitudes in Sweden. Arrows show the approxim- 
ate shift from univoltine to bivoltine phenology, a) Average for both generations in the bivoltine area, b) only the 
over-wintering (first) generation in all areas, c) only the directly developing (second) generation in the bivoltine area. 
Bars show standard errors. 



second generation. There were some individuals collec- 
ted between the two peaks that may belong to either of 
the tails of the two overlapping generations or in some 



cases represent a single generation in a cold year. These 
butterflies were excluded from relevant analyses. Since 
the generations differed in size in bivoltine areas, we 
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show results of depicting latitudinal patterns with these 
more southern areas represented by either all individuals 
(Figure 3a), only the first generation from over-wintering 
larvae (Figure 3b) or only the second, directly developing, 
generation (Figure 3c). In the field collection female size 
was consistendy smaller (Figure 3a— c; GLM ANOVA, n 
= 522, p<0.001); for this reason the sexes are shown 
separately. 

When all individuals from both generations were in- 
cluded there was no clear size pattern according to latit- 
ude (Figure 3 a). When only insects with a history of hi- 
bernation were included there was a pattern faindy re- 
miniscent of the predicted saw-tooth pattern (Figure 1), 
albeit with more smooth peaks and troughs (Figure 3b), 
and there was a significant effect of latitude on size 
(GLM ANOVA, n = 275 individuals, 10 latitudes, 
p<0.001; sex included as a factor). The butterflies were 
fairly large in the southernmost areas of Sweden, with 
size decreasing to the north. As predicted from the saw- 
tooth model, sizes increased again to a local maximum at 
around latitude 60° N, where the predominandy univolt- 
ine area begins, followed by a new decrease further 
north. Finally, size again increased at the very highest lat- 
itudes (Figures 3a— c). When only the second generation 
was chosen to represent the bivoltine areas (Figure 3c) 
the pattern was similar, but the maximum at latitude 60° 
N was even more evident. This is because wing length 
was always smaller for the direcdy developing generation 
than for the over-wintering generation in the bivoltine 
area below 60° N (compare Figures 3b-c; GLM 
ANOVA, n = 303 individuals, 4 latitudes, p<0.001; sex 
included as a factor). 

Common-garden experiment 

Development time did not differ between butterflies 
reared from individuals collected from the two low latit- 
ude sampling areas, i.e. those on the island of Oland (t- 
test, p = 0.84). Pupal mass did differ significandy (t-test, 
p<0.05), but variation among latitudes was much larger. 
For simplicity of analysis and illustration, these two 
southern sites were therefore pooled in all statistical ana- 
lyses and illustrations. 

Female butterflies reared from individuals collected from 
southern and intermediate areas had lower pupal mass 
than males, but this was not true for the strain from the 
high latitudinal areas (Figure 4a); this interaction was 
however not statistically significant (GLM ANOVA, two- 
way interaction between sex and latitude as categorical 
predictors, n = 304, p=0.095). When the interaction 
term was removed, the effect of sex was significant 
(p<0.001). If latitude was instead treated as a continuous 
predictor, the interaction was significant (p<0.05), as well 
as the effects of sex (p<0.001) and latitude (p<0.001). 
The heaviest pupae were produced in the southern strain 
and the lightest in the northern strain (Figure 4a). 



Interestingly, the latitudinal variation in mass did not 
correspond to variation in larval development time in the 
simple way assumed by Roff s theory on latitudinal pat- 
terns. Time spent in the larval stage differed significandy 
according to the latitudinal origin of strain (Figure 4b; 
GLM ANOVA, pO.OOl). However, larval development 
time was in fact shortest in the southern strain and 
longest in the northern, with the intermediate strain 
showing intermediate values (Figure 4b). There was no 
significant effect of sex on larval development time (p = 
0.96). Pupal development time differed significantly ac- 
cording to latitudinal origin (GLM ANOVA, p<0.05) but 
showing a pattern opposite to that for larval time (not 
shown). However, variation in pupal time was much 
smaller, so that the pattern for total development time 
was similar to that for larval time (effect of latitude 
p<0.001, not shown). The pattern of shorter develop- 
ment time yet larger mass in the southern areas corres- 
ponds to faster growth in southern strain, as illustrated in 
Figure 4c (GLM ANOVA, effect of latitude pO.OOl). 
The interaction between sex and latitude can also be 
seen for growth rate, but it is not significant (p = 0.46). 
Effect of sex was not significant with the interaction in- 
cluded in the model (p = 0.76), but when the interaction 
was removed a significant difference could be seen, with 
males growing faster than females (p<0.01). If latitude 
was instead treated as a continuous predictor, the inter- 
action was significant (p<0.05), as well as the effects of 
sex (p<0.01) and latitude (p<0.001). 

There was a significant effect of family (p<0.001 in all 
cases) in the analyses described above, if family was ad- 
ded as a random factor nested in latitude, strongly sug- 
gesting a genetic component behind the observed vari- 
ation. However, other results did not change; this is an 
indication that the geographical patterns may be robust 
even though there was a low number of families from two 
of the populations. 

In the Fl generation no correlation between develop- 
ment time and mass among individuals could be found in 
the low (n = 253, p = 0.9) and intermediate (n = 41, p = 
0.17) latitudinal sample, but in the high latitude strain 
they were positively correlated (n = 15, r^ = 0.36, 
p<0.05). 

Discussion 

The intention was to test if the saw-tooth pattern in size 
according to latitude and generation number, predicted 
by Roff (1980, 1983), is present in the common blue but- 
terfly P. icarus despite a lack of clear trade-offs between 
time and size in laboratory studies of this species (Leimar 
1 996). This is of particular interest in view of theory on 
three-dimensional time-rate-size relationships when 
growth rates are adaptively variable (Abrams et al. 1 996; 
Nylin and Gotthard 1998). Acommon-garden 
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Figure 4. Results of the common-garden experiment with Polyommatus icarus originating from three different latit- 
udes: Oland (Low=56.5° N), Stockholm (lntermediate=59.5° N) and Borgsjo (High=62.5° N). Arrows show approx- 
imate shift to a univoltine phenology, a) Mass (pupal weight), b) Larval development time, c) Growth rate calculated 
from the above. Bars show standard errors. 



experiment was also performed to obtain information on 
the genetic and plastic background of the phenotypic 
field pattern. 



Museum collection 

The geographical position of the shift from a univoltine 
to a bivoltine life cycle by P. icarus in Sweden is somewhat 
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ambiguous (see also Asher et al. 200 1 regarding a similar 
situation for populations in Britain). This can, at least 
partly, be explained by gene flow, preventing exact local 
adaptation. Climatic variation among years is another 
likely cause; early instar larvae will respond in a plastic 
manner to cues signalling how much time is left in the 
season, with different outcomes in different years. In 
warm years this will result in two generations and in cold 
years in a single generation, so that the "transition zone" 
(as seen in the combined museum material) will extend 
over large geographical areas. This is not simply an arte- 
fact of pooling years, but in a real sense represents the 
temporally fluctuating environment to which local popu- 
lations must adapt. In other words, most genetic differ- 
ences between areas must reflect the probability of a 
second generation in a given area, rather than a fixed life 
cycle pattern. 

The size of the P. icarus decreased towards higher latit- 
udes with shorter seasons in the bivoltine area, but then 
increased when they shifted to a univoltine lifestyle, ac- 
cording to predictions. Interestingly, this pattern was not 
evident when the generations were pooled in the bivolt- 
ine area (as is the usual procedure when phenology is not 
clear, e.g. Nylin and Svard 1991), because the genera- 
tions differed in size (see also Fric et al. 2006). This size 
difference in itself strongly suggests that adult size is, in 
part, determined by plasticity, since the major develop- 
mental pathway - diapause or direct development - fol- 
lowed by an individual is an outcome of high-level plasti- 
city according to seasonal cues such as photoperiod 
(Nylin and Gotthard 1998). 

Observations of this kind also suggest that size trends in 
insects may have been underestimated in the scientific lit- 
erature because of the difficulties introduced by genera- 
tion differences, and furthermore introduce a problem 
regarding which pattern should be used to test Roff s hy- 
potheses (1980, 1983). The pattern seen with only the 
first generation of adults (which has experienced larval 
hibernation) representing the bivoltine area (< 60° in Fig- 
ure 3b) could be argued to be the most correct test. This 
is because the univoltine individuals further north (> 60°) 
always over-winter and they are thus more comparable 
to the first generation of adults in the south. Hibernation 
is associated with special demands on an insect and could 
incur metabolic costs visible as reductions in adult size or, 
conversely, larger size if winter survival is positively size- 
dependent (Kemp and Jones 2001). 

On the other hand, it is the last generation of the favour- 
able season that has the most direct information on how 
much time remains before winter in the bivoltine areas, 
and if size is partly plastically determined this generation 
may show the clearest size responses as a result of end-of- 
season cues from the environment (Roff 1983). Plasticity 
theory (Leimar et al. 2006) suggests that in such 



ambiguous cases, information from the environment 
should affect what phenotype will result from an indi- 
vidual's genotype. In a sense, the genotype of an indi- 
vidual organism provides it with information on what 
was the fittest life cycle in past environments, but only in 
a probabilistic way, and in some situations the current 
environment provides more precise information that 
should be exploited if possible. 

Which generation should for these purposes be con- 
sidered the "last of the year", in the case of a bivoltine in- 
sect with late larval winter diapause such as P. icarus in 
Southern Sweden? This is not entirely clear either, but 
we suggest it is the second adult peak (Figure 3c). The 
first peak consists of individuals that spent the very last 
part of the previous summer as half-grown larvae and 
thus had a growth period at a later date than is ever the 
case for the second peak. However, the second peak con- 
sists of individuals that pupated later in the summer. 
Since growth of insects typically follows exponential or 
power functions, large size differences can occur due to 
small differences in the number of days spent in the last 
larval stage (Nylin 1994). This suggests that plastic size 
variation should mostly accumulate just before pupation, 
by the "decision" to continue growing or not as suggested 
by Leimar (1996) for P. icarus. The same conclusion can 
be reached based on the physiological model of body size 
variation in Lepidoptera presented by Davidowitz and 
Nijhout (2004), as they suggest that body size variation is 
the product of growth rate variation, coupled with vari- 
ation in critical mass for pupation, and in the interval to 
cessation of growth, once the critical mass has been 
reached. The fact that the second generation is consist- 
ently smaller in size supports this interpretation, since it 
may suggest that this generation has to "pay the cost" of 
a bivoltine pathway by pupating early to ensure that 
enough time in the season remains for its offspring to 
reach the hibernating stage. 

An unexpected pattern was found within the univoltine 
area. As predicted, sizes decreased going north within 
each of the areas of voltinism. However, further north in 
the univoltine area, P. icarus again increased its size, cor- 
responding to field observations (K. Fiedler, University of 
Vienna, personal communication) . If this pattern is to be 
explained by the saw-tooth theory, it would suggest that 
above latitude 65° P. icarus has adopted a two-year life- 
cycle Johansson 2003). However, a clear possibility re- 
mains that environmental factors and biological re- 
sponses other than those included in this theory come in- 
to play at such extreme latitudes. 

Common-garden experiment 

Regarding the laboratory experiments, note that the low 
number of individuals used to create the strains for the 
intermediate and high latitude populations means that 
we cannot claim to have a good picture of the genetic 
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make-up of these populations; more extensive sampling 
will be needed. Nevertheless, we believe that the differ- 
ences between geographical localities may prove robust, 
for two reasons. First, there was little difference between 
two sites at the more extensively sampled low latitude 
locality. Second, adding family as a factor to statistical in- 
vestigations did not change the overall results. 

The highly significant effects of latitudinal origin of strain 
demonstrate that there is a strong genetic component to 
the observed life history variation within Sweden, not 
only in size/mass but also in development time and 
growth rate. The potential for local adaptation is also 
suggested by the apparent presence of genetic variation 
for these life history traits, visible as significant variation 
among families. The presence of genetic latitudinal pat- 
terns in a common garden environment lead us to pro- 
pose that the life history observed in a given strain re- 
flects adaptation to the probability that the source popu- 
lation (averaged over years) will have enough time to 
complete two generations. This probability may well in- 
crease to lower latitudes in a near-linear fashion over the 
area covered by the three populations, and we propose 
that the decreasing development times and increased 
growth rates in the lower latitudes (Figure 4b— c) reflect 
this higher probability of a second generation. At the 
high latitude site a second generation is very rare or non- 
existent, hence these butterflies probably always have 
more than enough time to complete the single generation 
per year (at even higher latitudes the situation might be 
different, see below). Accordingly, it is likely that at this 
site the butterflies have been selected for slow growth, 
since fast growth could be cosdy (Gotthard et al. 1994, 
Gotthard 2000). Conversely, at the site from the lowest 
latitude it is likely that P. warns has been selected for fast 
growth, because there is generally enough time to com- 
plete two generations, but only barely. 

The strain developed from individuals collected at the in- 
termediate latitudes, were situated in a zone where devel- 
opment is not predictably uni- or bivoltine, and showed 
intermediate growth rates. This may be a result of selec- 
tion from an intermediate probability of a second genera- 
tion; however, gene flow from areas with more predict- 
able life cycles probably also contributes. The most ad- 
aptive life history in this situation could be very fast 
growth under direct development (not the observed inter- 
mediate rates, because they have even less time than the 
southern population to complete two generations) and 
very slow growth under development to diapause. Al- 
though differences between developmental pathways in 
this general direction are commonly found in butterflies 
and other insects (Nylin and Gotthard 1998) gene flow 
and other constraints on reaching the optimal reaction 
norm (such as genetic correlations among developmental 
pathways; Via and Lande 1985) may cause the 



intermediate growth rate and development time that we 
observed. 

It is interesting to note that although patterns of pupal 
mass according to latitudinal origin in the common- 
garden experiment followed a geographical pattern, the 
highest mass occurred at the lowest latitudes (Figure 4a). 
The field data show no such size pattern in either the first 
or the second generation, for these latitudes (compare 
Figure 3b-c with Figure 4a). Moreover, this pattern of 
higher mass at lower latitudes is not reflected by longer 
development times, the basic trade-off assumed by Roff 
(1980, 1983), but rather the opposite. Instead the more 
time-stressed butterflies at lower latitude grew faster 
(Figure 4c), achieving high mass in a short time. This can 
be seen as an example of counter-gradient variation 
(Conover and Schultz 1995), but not (as is typical) in rela- 
tion to field temperatures in the source populations, as 
this would have lead to the opposite growth rate cline. 
This is noteworthy because counter-gradient variation in 
response to temperature and time stress often co-vary 
and cannot readily be decoupled; this is the situation in, 
for instance, the common frog Rana temporaria where frogs 
from Northern Sweden grow faster (Laugen et al. 2003). 
It should be noted that the three sites for the common- 
garden experiment did not include such a time-stressed 
high latitude population, where there is just enough time 
for a single generation. At such sites we would expect 
laboratory growth rates to increase again, as in Alaskan 
Papilio canadensis butterflies when compared to Michigan 
populations at the same temperature (Ayres and Scriber 
1994). 

Combining the evidence 

We suggest that the discrepancy between field and labor- 
atory patterns in P. icarus can be ultimately explained by 
the fact that body size is the complex product of growth 
rate and development time (Wiklund et al. 1991; Abrams 
et al. 1996; Nylin and Gotthard 1998; Davidowitz and 
Nijhout 2004) and proximately by the constant laborat- 
ory environment. In the field, seasonal cues are available 
to signal the progress of the season and sometimes the 
need to cut development short at the expense of final 
size. In the common-garden experiment such cues were 
absent (other than a constant long photoperiod and high 
temperature) suggesting that such cues of high summer 
induced direct development in all individuals. Pupal mass 
was then free to reach some value that was, at least in 
part, a function of larval growth rates. Shorter photoperi- 
ods were not included within the scope of the present 
study, because the expected variation in diapause induc- 
tion within and among sites - with concurrent life history 
variation as noted above - would have added further to 
the complexity of interpretation. Short photoperiods can 
be expected to speed up growth in individuals destined 
for direct development, but not do so to the same extent 
in those destined for diapause (Wiklund et al. 1991; Nylin 
1992;Leimar 1996). 
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Similarly, the observed interaction between sex and site 
of origin for pupal mass (not seen in the field material) 
and growth rate demonstrates how observations of sexual 
size dimorphism in the laboratory may in part be the res- 
ult of variation in growth rate, rather than reflecting se- 
lection on mass per se. This suggests that the lack of sexual 
differences in growth rate and mass - or even reversed 
sexual difference - in the rearing of the northern popula- 
tion (Figure 4a and c) reflects release from selection for 
protandry under direct development in this univoltine 
population (Wiklund et al. 1991). Here, males can 
achieve protandry by simply breaking diapause early; 
hence in obligate univoltine populations there is no need 
for high larval growth rate in males. It has previously 
been shown for P. icarus (Fiedler and Holldobler 1992; 
Leimar 1 996) that the degree of sexual size dimorphism 
in this species is plastic. It has been demonstrated here 
that such laboratory results may not translate directly in- 
to size variation under natural conditions, although un- 
doubtedly they demonstrate the presence of sexual differ- 
ences in reaction norms for life history traits, and can il- 
luminate how life history decisions are taken according to 
sex. 

Despite evidence of growth rate variation in the laborat- 
ory, offsetting the correlation between final size and de- 
velopment time assumed by Rolf (1980, 1983), in the 
field we could still see a saw-tooth pattern fully consistent 
with his predictions. P. icarus is a relatively large insect re- 
stricted by season length and, hence, is predicted to fol- 
low the converse to Bergmann's rule, and its extension, 
the saw-tooth theory (Chown and Gaston 1999; Blanck- 
enhorn and Demont 2004). Evidently it does so in the 
field albeit not in a common environment. The grasshop- 
per Chorthippus brunneus provides an interesting comparis- 
on. In this species, time-constrained northern populations 
grow faster but have also evolved a lower critical mass 
(fewer instars) and thus shorter development time and 
lower adult mass (Telfer and Hassall 1999), in line with 
both Rolfs mechanistic assumptions and with his theor- 
etical predictions. In another grasshopper, Omocestus virid- 
ulus, populations from different altitudes grow at the 
same rate but again insects from time-stressed popula- 
tions have both a shorter development time, fewer instars 
and lower final mass in the laboratory (Berner and 
Blanckenhorn 2006). Future studies may reveal whether 
this stronger role for genetic variation in development 
time, as measured in the laboratory, is a general pattern 
explaining size variation in grasshoppers. An indication 
that this may perhaps not be the case is the fact that a re- 
view of intraspecilic variation in the number of larval in- 
stars in insects (Esperk et al. 2007) found that such vari- 
ation only rarely functions to couple development time 
positively with adult mass, as envisaged by Rolf (1983) 
and as found in these grasshoppers. 



Conclusions 

In conclusion, it is clear that the study of latitudinal and 
altitudinal patterns in size and other life history traits is 
far from complete. There is little agreement even on 
what is the most common pattern in a given group of an- 
imals, let alone on how the patterns are created by genet- 
ic and/ or plastic variation. As the investigation described 
here shows, the developmental processes behind the ob- 
served phenotypic patterns, even when they correspond 
to predictions from basic life history theory, can be more 
complex than assumed by such theory. Size differences 
between developmental pathways illustrate the role of 
phenotypic plasticity in these processes, but simultan- 
eously create problems for researchers already at the 
stage of establishing what the field and laboratory pat- 
terns really look like. To understand local life history ad- 
aptation we need both a larger database and more stud- 
ies where field data is combined with the results of 
common-garden experiments. 
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